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Abstract 

Gluon fusion into a very heavy neutrino pair by Higgs exchange is shown to 
lead to substantial production cross sections at pp supercolliders even without any 
extra generation of quarks. Rates are calculated for scalar as well as pseudoscalar 
Higgs. The angular correlation between dileptons emerging from the decays of 
the neutrinos shows distinctive features for Dirac and Majorana neutrinos as well 
as for scalar and pseudoscalar Higgs. 
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The production and subsequent decays of very heavy neutrinos (generically denoted 
with Mtv > Mz/'i) in forthcoming colhders has attracted a certain amount of recent interest 
[1-6]. Such particles are predicted by many theories — e.g. a 4th generation extension P] 
of the Standard Model (SM), left-right symmetric theories as well as 5*0(10) or i^g-based 
grand unified models P]. However, their properties are largely unconstrained by low energy 
phenomenology or astrophysical/cosmological considerations so long as they are unstable, 
stable ones being cosmologically disfavoured 0. Thus the best hope of detecting them rests 
on their copious production at high energy supercolliders and on the observation of their 
decays into distinctive final states. 

The requirement of a reasonable signal to background ratio in the above process points 
towards reactions induced by neutral currents rather than those effected by charged currents 
Taking this cue, we concentrate on the gluon fusion mechanism in a pp supercollider 
producing a Z boson or a spin zero Higgs particle (off-shell, in general,) which later goes 
into a pair of very heavy neutrinos that decay within the detector. The naive expectation 
would be for the quark-antiquark annihilation mechanism of Drell and Yan to produce the 
dominant contribution. Yet, this is belied |]^ on account of two reasons. (1) The density of 
gluons in the proton is far higher than that of quarks or antiquarks in the kinematic region 
of interest. (2) The quark mass appears as a factor in the coupling of quarks to the Higgs 
as well as to the longitudinal component of the Z suppressing the light quark contributions. 
Thus gluon fusion (through a heavy quark triangle) into a Z or Higgs and then into an A^A^ 
pair is the dominant production process at high energies. Ref has already discussed the 
Z-mediated process. Our focus in this communication is on the Higgs contribution to the 
gluon-fusion mechanism and some of its interesting properties. 

We assume that the heavy neutrinos decay within the detector into a light charged lepton 
(e, /i, r) by charged current interactions, i.e. N{N) — ^ + a real or virtual W converting 
into a lepton-neutrino pair or into jet(s). Universality constraints on such couplings are 



rather weak at present |12]. For heavy neutrinos belonging to a fourth generation, a general- 
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ized GIM mechanism would operate and these would be the only decay modes. For the case 
when Nl is an SU (2) x f/(l) singlet (and/or is part of an SU{2) doublet) there could also 
be a [0] neutral current induced decays: N{N) — > vl^l) + a real or virtual Z. However, 
each such decay would generate a very largely associated with a possibly reconstructible Z. 
Thus the corresponding events would have rather spectacular signatures. Our concern here 
is with the less spectacular but nevertheless interesting signals generated from the charged 
current induced decays of A^, A^. 

Thus, when is Majorana, we would have a like-sign dilepton pair without I^t and with 
jets, arising from the semi-hadronic decays of the pair. For a Dirac pair, though, one has 
to be more careful because of the background from the production and semi-leptonic decays 
of a tt pair. The best signal to trigger on now is the one consisting of the leptonic decays 
of both the heavy neutrinos. Despite the loss in terms of the branching fraction, one gains 
enormously in the cleanliness and uniqueness of the signal. The final state is characterized by 
four charged leptons, I^t and no jets. The choice of this hadronically quiet mode then, at a 
single stroke, gets rid of the potentially troublesome background from a ti production. With 
this constraint and the further realization that the 4-lepton signal from ti is suppressed by 
the branching fractions and that two of these leptons have, on the average, much lower E^, 
one sees easily that such backgrounds need not be worried about. The prominent background 
then arises from the production and decays of double Z's. With a cut on I^t and additionally 
a check on the unlike sign dilepton invariant mass distributions, this can be tackled. 

We assume a phenomenologically interesting mass-range of 50 - 500 GeV for a generic 
A^. The mass-range of the exchanged Higgs is also taken to be similar. For simplicity, 
consider either a pure scalar or a pure pseudoscalar coupling of the latter to the heavy 
neutrino pair. Such is the case in several interesting models — in particular — in the 
minimally extended SM with two Higgs doublets (with or without supersymmetry) . Left- 
right symmetric models, also have this property in consequence of the required absence of 
any tree level flavour-changing neutral current. We need not, however, confine ourselves to 
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any particular model at this stage and will come to specific mo del- dependence later. 
The above assumptions restrict the relevant Yukawa Lagrangian to 

Cy = {GfIV2)"^ J2 {{hs)iQrQiH + {hp),Qa5Q^A} + {h'sNNH + h'^N^.NA} 

Li 

(1) 

The index i here runs over the heavy quarks Q = {t, b). Terms involving the lighter quarks 
and ordinary leptons have been suppressed. H and A are scalar and pseudoscalar fields 
respectively while hs^p{h'g p) are real constants, their departure from unity refiecting any 
possible deviation from the SM couplings caused by either fermion mixing or Higgs mixing. 

In many models ( e.g. , one with two Higgs doublets), there are more than one neutral 
Higgs. Of these some have scalar couplings with the fermions, and we shall generically call 
them H. There are also those with pseudoscalar couplings and are generically designated A. 
Their vertices with two gluons, at one-loop order, is generated through heavy quark triangle 
diagrams. The gluon-gluon-scalar [0 vertex (with on-shell gluons of four momenta ki and 
k2, ki + k2 = q, but generally off-shell Higgs) is : 



S;t{ki, k2, q = k, + k2) = 6^' { k,,k2, - \g^, ] Y.{hs)m^ U 1 - - ) 



(2) 

In (0) = /m^ and} f{a) = /q^ dx ln[l — x(l — x) a + it\. Furthermore, as{q^) is the 
running QCD fine structure constant. Similarly, for the gluon-gluon-pseudoscalar vertex, we 
have 



Vltik,, fca; q = ki + k2)=i k^k^e^^^^ J2i^p)i m,/(a,). (3) 

T^q i 

A comparison of the above expression with that for the ggZ* vertex ^, ^ is instructive. 
The latter has an enhancement p] at high energies because of Yang's theorem, yet it is 
unable to dominate over the former. The key reason for this fact lies in the group theoretic 
nature of the QiQiZ coupling. The leading term in the ggZ* vertex comes proportional 
to TrlTsi) and hence vanishes as a part of the anomaly cancellation. The next-to-leading 
term is numerically much smaller. Such is definitely not the case for either of the ggH and 
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ggA vertices. Between the two, though, the pseudoscalar vertex is less affected by such 
cancellations in comparison with the scalar one. Thus, for equal coupling strengths, the 
pseudoscalar amplitude is somewhat larger. 

The differential cross section for the process g{ki) g{k2) — >N{pi, Si) N{p2, S2), with S'^ 

H*{q) 

standing for the spin pseudovector, is 
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In contrast, for giki) g(k2) — >N(pi,Si) N(p2,S2), it is 
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In (^) and denotes the subprocess cross section for producing a Dirac pair A^, out 

of gluon-gluon fusion. For a Majorana pair with similar couplings, one would instead have 
g] 0"*^ = 2(3"^. The observable cross-section a^pp NN + . . .) in the laboratory has to be 
obtained by folding the above expressions with the gluon distribution functions for which we 



use the GGR densities, as parametrized in Ref. |Tl| using the input of DFLM . 

The spin-summed integrated cross sections for the laboratory process pp NN + X have 
been plotted in Figs. (1 a,b), for various values of tuh- We have taken two representative 
cases. In Case I we have just three quark generations and have included only the t— and b— 
contributions with nit chosen to be IGOGeV (a change of 20GeV in nit would affect the cross- 
section at the level of only a few percent, an uncertainty inherent in the parametrizations 
for parton structure functions). In contrast. Case II includes an additional heavier fourth 
generation with nearly degnerate quark masses mu ^ 400 GeV. In order to be definite, 

we concentrate here on very heavy Dirac neutrinos and take the mixing factor at the HNN 
vertex to be unity. [This is the case not only for a fourth generation N, but also occurs for 
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those Eq GUTs where the same Higgs couples to both the SU{2)l x f/(l)y singlet neutrinos 
and the extra Q = —1/3 quarks present in the fermion 27-plet.] The corresponding plots 
for the pseudoscalar case occur in Figs. (2 a, b). The result could be viewed from a different 
perspective if the cross-sections were plotted as functions of mn instead (Fig. 3). This 
provides us with a ready estimate of the range of rriH for which a heavy neutrino of a given 
mass and coupling would be observable at the SSC. 

We want to emphasize the point that the cross sections are quite large even without 
the introduction of 4th generation quarks. In fact, for a very wide region of the parameter 
space, these are orders of magnitude larger than that for the Z-mediated process. This is 
true even when one is significantly away from the peaks (corresponding to the s-channel 
Higgs resonance) in Fig. 1-3. This result, for Majorana neutrinos (since only those were 
considered there), was basically contained in the work of Datta and Pilaftsis 0. However, in 
all realistic seesaw models of Majorana mass-generation, there is a small mixing angle in the 
HNN coupling. This led to an effective suppression of the cross section in their case. Thus, 
both in Refs. [|5|, |^, an extra generation of heavy quarks had to be brought in to generate 
measurable rates. Such a suppression factor need not be present for a heavy Dirac neutrino. 

An additional feature of our investigation is the pseudoscalar channel. Given the not 
unnatural assumption of comparable heavy neutrino coupling strengths for the scalar and 
the pseudoscalar Higgs, it is evident from Fig. 1-2 that the cross section for the pseu- 
doscalar exceeds that for the scalar by almost an order of magnitude. This is a new result 
of considerable significance — especially in the context of two Higgs doublet models (with 
or without supersymmetry) . Of course, this observation gets vitiated if the Higgs has simul- 
taneous scalar and pseudoscalar couplings; but the latter is incompatible with the principle 
of natural flavour-conservation [|l^ and we ignore it. 

The next issue, once the measureability of the heavy neutrino pair-production rate has 
been established, is the signal profile. As explained earlier, Majorana A^A^ production has 
the very distinctive signal of a pair of like-sign dileptons and jets but no I^t ■ In contrast, 
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Dirac A^A^ production is best studied through two pairs (one hard, one softer) of unhke- 
sign dileptons and nonzero I^t ■ The Et- and the rapidity-distributions of the charged 
leptons obviously follow two patterns: one for the lepton directly produced at the A^-decay 
(primary) vertex and the other from W^-decay. While the former distribution depends on 
the mass difference M^r — M^, the latter is much less dependent on and peaks around 
Et ~ 60 GeV. 

The double Z background can obviously be bothersome in the light of significant uncer- 
tainties possible in i?T-measurements at supercolliders. But, an acceptance cut of Et > 30 
GeV would eliminate most of it without affecting the signal much. However, if rriH < 150 
GeV, this last statement does not hold and one can lose a large part of the signal in impos- 
ing an ET-cut. A better method of reducing the background might be to reconstruct the 
Z-events. A study of the invariant mass distributions of a pair of oppositely charged final 
state leptons in the signal shows that these are not peaked around m£+£- ~ Mz- Thus the 
imposition of an acceptance cut of |m^+£- — Mz \ ~ 5Tz would reduce the signal strength only 
by about 15 — 20% for mi^ 100 — 200 GeV . On the other hand, this disposes of almost 
the entire ZZ — > background. 

Once the heavy neutrino has been detected, the most obvious parameter to be determined 
next is its mass. This can be done most effectively by looking at the m£+£- distribution in 
each hemisphere. For a sufficiently massive A^ decaying into three massless leptons through 
charged current interactions, the VT-resonance would dominate by far. For such a case, the 
above invariant mass distribution would have a very sharp cut-off at m^+i- = (m^ — m^)^/^, 
thus determining rriN- 

Having established the feasibility of searching for such heavy neutrinos, let us now turn 
to the proposal of distinguishing between Dirac and Majorana ones, first mooted in Ref. 0. 
For the large values of M^, considered here, the measurement of final state lepton charge 
signs may not be easy since these leptons are going to be very energetic. Nevertheless, the 
angular correlation between two such leptons should be measurable. The authors of Ref. [Q 
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employed the formalism of Tsai [|T5| to show that the dilepton angular correlation for the 



process gg—^NN, N iX, N — > J'X' {X, X' unobserved) discriminates between Dirac and 
Majorana A^'s. Such a procedure works in the present case too. This is straight forward to 
see once one realizes that the discrimination in Ref. |P in fact originates from the contraction 
of the ggZ* vertex and the Z*-propagator to the ZNN coupling which effectively renders the 
latter into a pseudoscalar one.^ It follows that, for a pseudoscalar exchange, the correlation 
is identical to that of Ref. . For the scalar-mediated case the structure of the vertex is a 
little different, but nonetheless the discrimination works. 

In order to elaborate on the last statement, we go to the gg CM frame and write the 
momenta and spin pseudovectors of and (with velocities ±p) respectively as 

Vl,2 = {E, 0, 0, ±(3E), S,,2 = [±I3{1 - ei,2 6,2 y, (l ' . 

where are the spin-vectors in their corresponding rest frames. The relations 



2), 



{l-S,-S2){pi-p2-Ml)+p,-S2P2-Sr = 2E^(3\l+ir-i2) 
follow. The use of these relations enables us to rewrite (^) and 
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Mq') (1 + ^1-6), 

-V 
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dVlj^f 

where we have lumped all the ^independent parts into the factored function fH,A{<f)- 

We can further designate the polarization vectors^ of and A^ as uji and UJ2 respectively 
and denote the momenta of the final state decay leptons ^(?), for A^ iW, N iW) by 
qi{q2)- Then, in the notation of Ref. ||^, we can write 

dTjN ^ iX) ^ - -^ 

dT{N ^iX) ^ , ^ . . , 

— = CAr(a7v + /^/Arga ■ ^^2), 

d\ If 
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with 



B being the branching ratio for the W-dec&y that is being triggered on ( e.g. in the Dirac 
case it would be the leptonic branching fraction). 

Again, as in the use of Eq. (4.26) of Ref. [1^ leads using eqs. (0) and (|D to the 
triple differential cross sections for a Dirac 

,o wo = WD 12 ( "TV + I^nQi ■ (l2 ) and 

dilNdilidilj \}-tot\ ^ ' /.pN 

with r^^Q^" as the total decay-width of the Dirac A^. In the Majorana case if the lepton 
charge signs are unobserved, one must take the sum 

^3^M ^3^M ^3^M 

+ ,o • (11) 



dQjydflg-dQi- dQ]\fdQe+dQi+ dQ]\fdfli+dfl 
This can be calculated following Ref. |^ to be 



d'a''{H/A) 16fH,Ai 



^2> 



CWn- (12) 



dVLj^dVLji^dVtii^i ^^tot) 
In ([T2|) r|,Qy is the total decay width^ of the Majorana A^. 

The primary charged dileptons emanating from a pair of decaying heavy Dirac neutrinos 
do show a distinctive angular correlation depending on whether the process involves scalar 
or pseudoscalar exchange. For the latter this correlation is identical to that in the Z- 
exchange case. Exactly, as in that case, here too the charged leptons do not show any 
angular correlation when A^ is a Majorana particle. In reality, all three exchanges (Z, H and 
A) will contribute to the process gg A^A^. But these different contributions correspond to 
different initial state angular momenta and parity, i.e. different partial waves, and hence 
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add incoherently. Thus the basic dilepton angular correlation discrimination — yes for Dirac, 
no for Majorana — remains intact. There is the additional novel result in this work that 
the angular correlations for the Dirac case are different [vide eqns. (p!0D ] for the scalar and 
pseudoscalar exchange contributions. To be sure, these correlations have been discussed 
in the gg CM frame. However, given the gluon density distributions, the transformation 
of these results to the pp CM. frame in the laboratory can be done by straightforward 
Monte-Carlo methods. 

Though the above argument holds for N{N) — > £(I) + X, for all X, in practice, for 
hadronic states X, the background is very large in a pp collider. Thus one has to confine 
oneself to purely leptonic X. But this could lead to a possible contamination in the correla- 
tion as the argument works only for the charged leptons coming directly from N{N). One 
could of course get around this problem if the two leptons coming from the same heavy neu- 
trino could be distinguished. Once the mass of the very heavy neutrino has been determined, 
PT^distributions coupled with kinematic considerations then do this job for us. Of course 
since we are talking of three body decays, one can never completely isolate the secondary 
lepton, but for sufficiently heavy N, there is very little overlap between the pt distributions. 

In summary, we have presented a feasibility study for finding very heavy neutrinos at pp 
supercolliders through direct pair-production from gluon fusion. For the mass-range that is 
kinematically interesting at such machines, this process dominates. We have compared the 
contributions due to different propagating bosons in the intermediate state arriving at the 
following conclusion. For a large area of the parameter space, Higgs-exchange overwhelms 
the Z-contribution, with the pseudoscalar (if existing) taking the lead. This is to the extent 
that measurable supercollider rates are possible (because of the heavy top) with only three 
quark generations — a fact that is not true for the Z exchange contribution. We have 
also shown how final state dilepton angular correlation can discriminate a Majorana heavy 
neutrino from a Dirac one and, for the latter, a scalar Higgs exchange contribution from a 
pseudoscalar or Z-mediated one. Much of what we say also goes through for the production 
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of neutral fermions in a supersymmetric theory. These and other related issues will be 
elaborated elsewhere. 

We thank A. Datta, M. Drees and D.P. Roy for helpful discussions. 
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Footnotes 

1. The function /(a) = 2 (sin"^ ^faji) for < a < 4 and 2 (cosh"^ ^0/2)^ - vrV2 + 
2i7rcosh'^ -^/a/^ for «- > 4. 

2. In effect, (gi + g2)^[-fi'pA + (gi + g2)p(gi + q2)\Mz'^\ un{Pi)1i,{9v + fl-ATs) ^^7v(P2) = 
2MNgAMz\S - M|) M7v(pi)75yiv(P2). 



3. 'Wii{w2i) is JTSl the ratio of the difference between the number of A^'s (A^'s) polarized 
along the positive ith axis and that of A^'s (A^'s) polarized along the negative ith axis 
to the sum of these two numbers. 



4. If the field A^ couples only to the lepton field £ and the charged- weak-boson field W, 
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Figure Captions 

1. a{pp — > NN) mediated by a scalar Higgs and without mixing suppression. 

a) Only three generations of quarks with = 160 GeV 

b) Including an additional fourth generation with rrit' = my = 400 GeV. 

The continuous, dotted, dashed and the dot-dashed lines correspond to ttih — 100, 200, 300, and 600 C 
respectively. 

2. Same as Figs. la,b, but with a pseudoscalar Higgs exchange instead. 

3. Same as Fig. la but with mn as abscissa. 

The continuous, dotted and the dashed lines correspond to uin — 100, 200, and 300 GeV 
respectively. 
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